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I review recent work on the calculation of hadron masses, decay constants and wave functions. 



1. INTRODUCTION 

Exactly what set of topics is supposed to go 
into the spectroscopy review talk seems to vary a 
bit from one lattice conference to the next. My 
definition of the territory this year is masses, de- 
cay constants and wave functions for both the 
valence (quenched) approximation and full QCD. 
I will include glueballs but not hadrons contain- 
ing heavy quarks and will restrict myself mainly 
to work done during the last year. 

The spectroscopy of low lying quark states and 
glue states I think remains a critical part of the 
QCD. An unavoidable component of showing that 
QCD actually does explain the strong interac- 
tions is to show that it accounts for the masses 
of low lying hadrons. QCD has been assigned the 
homework problem of hadron spectroscopy just as 
early quantum mechanics was assigned to work on 
the spectrum of hydrogen and other light atoms. 
In addition, whatever algorithms can prove them- 
selves in reproducing the known features of low ly- 
ing spectroscopy can then be applied with greater 
confidence to extract predictions from QCD for 
things we do not know yet. The masses of the 
lightest glueballs, of course, are not known. A re- 
liable lattice calculation of these numbers seems 
likely to be a required first step toward the ex- 
perimental identification of glueballs. 

The overall picture of the present state of 
hadron spectroscopy is roughly as follows. For 
spectroscopy including quark-antiquark vacuum 
polarization, life is still pretty difficult. Calcu- 
lations are restricted to fairly small statistical 
ensembles and comparatively large lattice spac- 
ing and quark mass. Significantly more machine 
power, perhaps a factor of 100 or more, or a 



real improvement in algorithms, are needed to 
get reliable estimates for the zero lattice spac- 
ing, infinite volume limits of most properties of 
light hadrons. For full QCD calculations and 
valence approximation calculations tuned to the 
same values of rn^^/rnp, lattice spacing in physi- 
cal units and lattice period in physical units, and 
with reasonably large lattice period, correspond- 
ing hadron masses, at the parameter values for 
which calculations have been done so far, agree 
within statistics. Full QCD and valence approx- 
imation calculations, with parameters tuned to 
agree at sufficiently large volume, yield masses 
which disagree if the the lattice volume is then 
made much smaller. The difference in volume de- 
pendence between full QCD and the valence ap- 
proximation has been shown to be a consequence 
of the nonzero vacuum expectation value which 
occurs in full QCD for gauge loops closed around 
a space-direction lattice period. Except for very 
small lattices, these loops have zero expectation 
in the valence approximation. 

For spectroscopy in the valence approximation, 
on the other hand, there are now calculations for 
Wilson quarks of masses and decay constants ex- 
trapolated to physical quark mass, zero lattice 
spacing and infinite volume. Ratios of several 
light hadron masses found in these limits dif- 
fer from experiment by less than 6%, with sta- 
tistical errors of 8% and less, and are statisti- 
cally consistent with experiment. Although chi- 
ral perturbation theory suggests that the valence 
approximation will behave anomalously at suf- 
ficiently small quark mass, it now appears that 
this problem would only disrupt the hadron mass 
ratios which have been calculated if they were 



evaluated at quark masses below the light quark 
masses. Valence approximation meson decay con- 
stants extrapolated to physical quark mass, zero 
lattice spacing and infinite volume fall typically 
about 15% below experiment. This underesti- 
mate ranges in significance from about 1 to 3 
standard deviations. The valence approximation 
is expected, however, to be most reliable for quan- 
tities determined primarily by the low momentum 
behavior of the chromoelectric field. Decay con- 
stants are proportional to the absolute square of 
meson wave functions at the origin and are thus 
sensitive to the high momentum behavior of the 
chromoelectric field. A simple renormalization 
group estimate suggests that the error in the va- 
lence approximation's treatment of the high mo- 
mentum behavior of the chromoelectric field will 
lead to decay constants falling below those of the 
full theory. 

Two independent, infinite volume, continuum 
limit valence approximation calculations have 
now also been reported for the scalar glueball 
mass. The two data sets agree within statistical 
errors before extrapolation, but the two groups 
extrapolate differently to zero lattice spacing and 
get answers which are not quite consistent. The 
calculation with higher statistical weight, using 
about eight times as many configurations, pre- 
dicts a mass of 1740 ± 70 MeV, favoring /o(1710) 
as the scalar glueball. 

The subjects I am going to cover will be orga- 
nized as follows: 

1. Technical issues 

• Extracting hadron masses from prop- 
agators 

• Volume dependence 

• Valence approximation versus full 
QCD 

• Wave functions, state vectors 

2. Calculations with fixed lattice spacing 

• Valence approximation 

• Full QCD 

3. Valence approximation zero lattice spacing 
limits 



• Masses 

• Decay constants 

• Glueball masses 

2. TECHNICAL ISSUES 

2.1. Extracting Masses from Propagators 

At Lattice 92 the QCDPAX collaboration [1] 
presented calculations checking for high mass con- 
tamination in the values of hadron masses re- 
ported by several other groups. The lightest mass 
contributing to a propagator should be extracted 
from the propagator's fall off at asymptotically 
large time separations, and picks up high mass 
contamination from excited intermediate states 
if it is obtained at time separations which are 
too small. The value which QCDPAX found 
for the rho mass with Wilson quarks in the va- 
lence approximation using 200 gauge configura- 
tions on a lattice 24^ x 54 at /3 of 6.0 and k of 
0.1550 was 3 standard deviations below the value 
which the APE collaboration found at the same 
/3 and k using 78 gauge configurations on a lat- 
tice 24^ X 32 [2]. Both groups have now run new 
calculations with larger ensembles [3,4]. The new 
calculations of the rho mass are both in agree- 
ment with the QCDPAX value from Lattice 92. 
The original disagreement appears to me to have 
been caused by a statistical fiuctuation in the first 
APE ensemble. The fiuctuation caused the ap- 
pearance of an effective mass plateau at a some- 
what smaller time separation and therefore with 
a somewhat larger mass value than occurs in the 
new data. 

At Lattice 92 Ukawa [5] questioned whether 
high mass contamination might not be present 
in some of the masses which I discussed from 
calculations on GFll [6]. A comparison of 
masses found from propagators for several dif- 
ferent choices of hadron sink operators, however, 
provides evidence that our numbers probably do 
not have significant high mass contamination. 
This subject will be discussed below in Sect. 4.1. 

2.2. Volume Dependence 

New data, extending earlier results [7], on the 
volume dependence of hadron masses in full QCD 
was presented by the MILC collaboration [8]. 
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Figure 1. Volume dependence of hadron masses 
found by the MILC collaboration at /3 of 5.415. 
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Figure 2. Volume dependence of hadron masses 
found by the Kyoto-Tsukuba collaboration at /3 
of 5.7. 



These calculations were done with full QCD for 
two flavors of Kogut-Susskind fermions. A sum- 
mary of the run parameters is shown in Table 1. 
For the range of parameters examined, values 
of TOjv and rUp within 2% of their infinite vol- 
ume limits required space direction lattice periods 
greater than about 2.5 fm. The volume depen- 
dence of hadron masses for rUqa of 0.05 is shown 
in Figure 1. The horizontal lines indicate inter- 
vals within 2% of infinite volume masses. The lat- 
tice spacing was chosen by setting the pion mass 
to 500 MeV, since in this case rriT^/mp is nearly 



0.597(1) 1.023(2) 
0.489(2) 0.918(4) 
0.367(3) 0.883(6) 



Table 1 

Parameters of MILC collaboration runs, on lat- 
tices S^ X 24 with rriT^/mp in each case given for 
the largest lattice. 



/3 


rUqa S mT^/mp 


rupa 


valence approximation 


6.0 
6.0 


0.02 6-24 0.643(4) 
0.01 6-24 0.513(7) 


0.520(3) 
0.465(6) 


full QCD 


5.7 
5.7 


0.02 8-20 0.692(5) 
0.01 8-20 0.586(11) 


0.492(3) 
0.418(7) 



Table 2 

Parameters of runs by the Kyoto-Tsukuba collab- 
oration, on lattices S^ x 24 with mT^jmp and rripa 
in each case given for the largest lattice. 

equal to the observed value of to^/to^. 

Data on the volume dependence of valence ap- 
proximation masses [9] and a comparison with the 
volume dependence of full QCD masses [10] was 
reported by the Kyoto-Tsukuba collaboration . 
The parameters of these runs are shown in Ta- 
ble 2. For the largest values of space direction 
period, the values of to.^/top and rripa for the va- 
lence approximation run at TOjO of 0.01 appear 
roughly comparable to those of the full QCD run 
at the same quark mass, and a similar approxi- 
mate equality holds for the runs at TOjO of 0.02. 
Data for full QCD and for the valence approxima- 
tion with TOjO of 0.01 is shown in Figure 2. In the 
valence approximation, however, rn.^, rrip and tojv 
are found to approach their infinite volume val- 
ues more rapidly than they do in full QCD. The 
origin of this difference is shown to be related to 
the difference between the behavior in full QCD 
and in the valence approximation of the variables 
Pi, the trace of the product of gauge links closed 
around the lattice period in space direction i. The 
values of Pi affect the contribution to quark prop- 
agators of quark paths wrapping around the lat- 
tice period. In the valence approximation, for 



lattice periods larger than the inverse of the de- 
confining temperature, the familiar Z^ symmetry 
of the QCD action without fermions prevents Pi 
from acquiring a vacuum expectation value. In 
full QCD this symmetry is broken by the fermion 
terms in the action, and Pi has a nonzero vacuum 
expectation value which falls continuously to zero 
as the lattice period becomes larger. By changing 
the boundary conditions of the fermions entering 
the action from periodic to antiperiodic, the sign 
of the vacuum expectation of Pi can be changed. 
The Kyoto-Tsukuba collaboration shows that this 
change causes hadron masses to go from decreas- 
ing with volume to increasing with volume. This 
result, along with similar information obtained 
varying the boundary conditions on the quarks 
entering hadron propagators, establishes that the 
difference between the vacuum expectation of Pi 
in the valence approximation and in the full the- 
ory accounts for the differences in volume depen- 
dence. 

The volume dependence in hadron masses in 
full QCD reported by the MILC collaboration and 
by the Kyoto-Tsukuba collaboration are consis- 
tent with each other. For all but the largest values 
of space direction lattice period, both data sets 
shown significantly stronger volume dependence 
than predicted by Luscher's rigorous asymptotic 
formula [11]. Both sets of results are consistent 
with the behavior 



'i{L) = to(oo) + 



L^ 



(1) 



predicted for intermediate values of L by the 
model of Ref. [12]. Luscher's formula embodies 
the effect of meson exchange between a hadron 
and its nearest periodic image. The model of 
Ref. [12] is intended to take into account in ad- 
dition the finite extension of a hadron's wave 
function, or equivalently the momentum depen- 
dent form factor appearing in a hadron's cou- 
pling to exchanged mesons. For lattice peri- 
ods significantly larger than the extension of a 
hadron's wave function, the model of Ref. [12] is 
designed to reproduce Luscher's formula. Fits of 
the Kyoto-Tsukuba data for full QCD to Eq. (1) 
are shown in Figure 2. 



2.3. Valence Approximation versus Full 
QCD 

The valence approximation may be viewed 
as replacing the momentum-dependent color di- 
electric constant arising from quark-antiquark 
vacuum polarization with its zero-momentum 
limit [13]. This approximation might be expected 
to be fairly reliable for low-lying baryon and me- 
son masses determined by the low momentum be- 
havior of the chromoelectric field. 

Missing from the valence approximation, how- 
ever, are couplings of vector mesons to their de- 
cay channels. For the spin-3/2 baryon multiplet, 
these couplings are present but altered from their 
values in full QCD [14]. A calculation by Leinwe- 
ber and Cohen [15] attempts to estimate the effect 
of this omission on irip. These authors consider 
also the effect of obtaining irip by the linear ex- 
trapolation down to light quark masses of values 
of vector masses calculated with heavier quarks. 
Leinweber and Cohen use a simple local coupling 
between the rho and two pions, combined with a 
momentum cutoff in the otherwise divergent pion 
loop integral entering the rho propagator. The 
cutoff is intended to model the momentum de- 
pendence of the coupling between the rho and pi- 
ons which, in the real world, would cause the loop 
integral to converge. The error in the valence ap- 
proximation value for irip arising from the missing 
loop combined with extrapolation from the range 
of quark masses used in Ref. [6], is estimated to be 
between and —25 MeV, thus less than 3%. For 
a linear extrapolation of the vector meson mass 
as a function of valence quark mass in full QCD, 
working from a quark mass interval extending up- 
ward from the point corresponding to a 500 MeV 
pion, the error arising from extrapolation alone is 
estimated to be -10 to -20 MeV. 

A discussion in Ref. [6] of the extrapolation of 
hadron masses produces as a byproduct an addi- 
tional piece of evidence suggesting that the cou- 
pling of light hadrons to decay channels has a rel- 
atively small effect on their masses in full QCD 
and therefore gives rise to a relatively small error 
when omitted or altered in the valence approx- 
imation. Take the up and down quark masses 
to be equal, for simplicity, to a single normal 
quark mass, to„. Then expand the masses of 




Figure 3. Linear extrapolations to determine 
nonstrange vector, spin-1/2 and spin-3/2 baryon 
masses from the masses of their strange partners. 

the hadrons in the vector, spin-1/2 and spin-3/2 
baryon multiplets as power series in the mass 
of the strange quark, rn^, around the point at 
which the strange quark mass and normal quark 
mass are equal, keeping only the leading linear 
term. Applying flavor SU(3) symmetry at the 
point with rUs equal to to„ then gives relations 
between the masses of hadrons composed of a sin- 
gle flavor of quark and hadrons composed of two 
different flavors of quarks. For the vector multi- 
plet we obtain, for example. 



mp[{mn + TOs)/2, (m„ + ms)/2] = 
mK'{mn,m,). 



(2) 



In addition, in the valence approximation, m^ is 
given by TOp(tos,TOs). Corresponding sets of re- 
lations hold in the baryon multiplets. 

Using these relations, the mass of the member 
of each multiplet containing no strange quarks 
can be obtained by linear extrapolation from the 
masses of the multiplet members containing one 
or more strange quark. Extrapolating from ex- 
perimentally observed masses in this way requires 
only the quark content of each hadron involved 
and does not depend on the actual values of 
quark masses. These extrapolations are, essen- 
tially, just a reinterpretation of the Gell-Mann- 
Okubo mass formulas. Figure 3 shows extrapo- 



lations for the vector mesons, spin-1/2 baryons 
and spin-3/2 baryons. Leinweber and Cohen's 
calculation for full QCD may be viewed as an es- 
timate of the error of the extrapolation for the 
vector multiplet. In place of their error value of 
— 10 to —20 MeV, the error actually obtained ap- 
plying this extrapolation to the mass data from 
experiment is —4 MeV, or about 0.5%. For the 
spin-1/2 baryon multiplet the error is +13 MeV, 
about 1.4%, and for the spin-3/2 baryon multi- 
plet the error is +10 MeV, about 0.8%. In two 
of these multiplets, however, decay widths vary 
drastically among the strange hadrons to which 
the linear fits are made and are generally much 
smaller than the decay width of the nonstrange 
hadron whose mass is found by the extrapolation. 
In the vector meson multiplet, the rho width is 
149 MeV, while the K* and phi used to find the 
rho mass have widths of 49.8 MeV and 4.4 MeV, 
respectively. The members of the spin-1/2 baryon 
multiplet are stable with respect to the strong in- 
teraction. In the spin-3/2 baryon multiplet, the 
A has a width of about 115 MeV, while S*, S* 
and Q used to find the A mass have widths of 
36 MeV, 10 MeV and MeV, respectively. The 
Q does not decay by the strong interaction. It is 
hard to imagine how the masses of the broad rho 
and A could be extrapolated so accurately from 
their strange partners, whose widths are much 
narrower and vary significantly, if hadron masses 
were altered by more than a percent or so by their 
decay widths. 

Calculations of the difference between the va- 
lence approximation and full QCD in the limit 
of small quark masses have been reported by 
Bernard and Golterman [16] and by Sharpe [17]. 
For hadrons composed of quarks which are suffi- 
ciently light, these authors show that chiral per- 
turbation theory suggests the valence approxima- 
tion will behave pathologically. As a consequence 
of the absence of vacuum polarization, the rj' mass 
goes to zero along with the pion mass as quark 
masses are made small. Virtual rj' loops become 
divergent at low momentum as the quark mass is 
made small and contribute additional logarithms 
of the quark mass not present in chiral perturba- 
tion theory for full QCD and not multiplied by 
as many powers of rn^ as would occur in the full 



theory. For the behavior of the pion mass at small 
quark mass, which in full QCD is 

ml = crug, (3) 

Sharpe has summed a collection of such singu- 
lar diagrams for the valence approximation and 
obtained 



(4) 
(5) 
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Here /j, is a measure of the amplitude for the 
quark-antiquark pair in the rj' to annihilate into 
gluons. In full QCD /j, becomes the mass of the 
rj' in the limit of zero quark mass and would be 
about 900 MeV. Using this estimate, along with 
93 MeV for /^, the value oi 6 becomes about 0.2. 
There is, however, no reason to assume /j, has 
the same value in full QCD and in the valence 
approximation. A recent Monte Carlo evaluation 
of /J, in the valence approximation [18] yielded 
0.700 ± 0.050 GeV, giving 6 of about 0.12. Fits of 
Eq. (4) to data to be discussed in Sect. 4.1 for 
Wilson quarks [6] and for Kogut-Susskind [25] 
quarks, give values of 6 shown in Table 3. For 
Wilson quarks, rriq in these fits is (2k)~^ — (2kc)~^ 
where k^ is the critical hopping constant at which 
TO^ becomes zero. The value of k^ was taken as 
one of the adjustable parameters in each fit. The 
negative values of 6 for two out of three of the fits 
with Wilson quarks imply that the corresponding 
rUq are above the range in which the argument 
leading to Eq. (4) is applicable. The values of 
6 found with Kogut-Susskind quarks are much 
smaller than 0.12 but consistently positive. The 
Kogut-Susskind fits are done at smaller rriT^/mp 
than those with Wilson fermions, and therefore 
would correspond to smaller values of rUq if mea- 
sured in common units. Thus there is some ev- 
idence for a slow approach to the behavior of 
Eq. (4), perhaps with 6 of 0.12, for rriT^/mp suf- 
ficiently below 0.3. In any case, if the extrap- 
olations to small quark mass to be discussed in 
Sect. 4.1 were done using Eq. (3) for rriT^/mp 
above 0.3 and using Eq. (4) with 6 of 0.12 for 
TO^/TOp below 0.3, in place of the actual fits us- 
ing Eq. (3) everywhere, the changes in predicted 



Wilson 



16^ X 24 

243 X 36 

30 X 322 X 40 



5.70 0.49 - 0.69 0.04(6) 
5.93 0.47 - 0.74 -0.03(3) 
6.17 0.48 - 0.74 -0.07(3) 



Kogut-Susskind 



6.00 
6.50 



0.31 
0.40 



0.50 
0.66 



f«0.03 
0.02(1) 



32^ X 64 

32^ X 64 
Table 3 

Values of the parameter 6 obtained from fits to 
data for the pion mass as a function of quark 
mass. 

hadron mass ratios would be smaller than the sta- 
tistical errors. 

Labrenz and Sharpe [19] have also attempted 
to use chiral perturbation theory to estimate the 
effect of the valence approximation on the nucleon 
mass, TOjv. While I believe that the basic idea of 
trying to use chiral perturbation theory to make 
such comparisons is very promising, the particu- 
lar results which they presented in Dallas appear 
to me to be unconvincing. For the behavior of 
TOjv at small pion mass, chiral perturbation the- 
ory for full QCD gives an expression of the form 



mjv 



b^ml +bKm%; +c^ml +CKm%;. (6) 



Chiral perturbation theory for the valence ap- 
proximation gives an expression with an addi- 
tional pion term and no K terms, since these en- 
tail vacuum strange quark loops. 



mjv 



bl 



(7) 



The extra a"TO^ term arises from the same mech- 
anism leading to Eq. (4), and the coefficient a" 
is proportional to 6. Our earlier discussion of 6 
suggests that the a" coefficient will be small. 

Although it might seem that a comparison of 
Eqs. (6) and (7) could be used to estimate the er- 
ror in valence approximation calculations of tojv, 
it turns out this can not be done. One difficulty is 
that Eq. (6) describes the variation of tojv as both 
sea and valence quark masses are varied together, 
while Eq. (7) describes only the result of vary- 
ing the valence quark mass with the sea quark 
mass kept fixed. As the sea quark mass varies 
in full QCD, the effective low-momentum gauge 



coupling varies, and therefore the gauge coupling 
to be put into the corresponding valence approx- 
imation must vary. Thus for each choice of rn^ 
and rriK in the full QCD relation Eq. (6), a dif- 
ferent and a priori unknown value of a' must be 
used in the corresponding valence approximation 
relation Eq. (7). For lattice spacing sufficiently 
small, a change in the gauge coupling entering 
the valence approximation causes a change in low 
lying hadron masses by a single overall scale fac- 
tor. Thus the ratio of Eq. (7) to a corresponding 
equation for some other mass parameter, such as 
/,r, will be free of the scale factor and perhaps 
might be compared to the ratio of Eqs. (6) to an 
expression for /^ in full QCD. 

A useful comparison for the ratio tojvZ/tt still 
can not be made. The difficulty at this point 
is that the valence approximation chiral pertur- 
bation theory expression for TOjv//7r depends on 
an unknown additive constant, comparable to a' 
in Eq. (7), giving the valence approximation to 
"^jv/Ztt in the limit of zero quark mass. I know 
of no way to fix this number from chiral pertur- 
bation theory. In the limit of very heavy rn^ and 
rriK, h^jv/Ztt will become equal in full QCD and 
in the valence approximation. This limit, how- 
ever, is beyond the range of applicability of chiral 
perturbation theory and provides no help in de- 
termining the missing parameter. One might try 
to determine the missing coefficient by a fit to the 
data of Ref. [6] . Unfortunately, as mentioned al- 
ready in the discussion of Eq. (4), this data does 
not appear to extend to sufficiently small quark 
mass for the full apparatus of chiral perturbation 
theory to be applicable. 

If, nonetheless, the expression of Eq. (7) is fit- 
ted [19] to the data of Ref. [6] at j3 of 5.93, the 
extrapolated value of tojv at the light quark mass 
differs from the result of the linear extrapolation 
used in Ref. [6] by only about 3%. 

There are comparisons between full QCD and 
the valence approximation which have fewer dif- 
ficulties. An example of one such comparison is 
for /kt/Ztt. In both full QCD and in the valence 
approximation this quantity must become 1 if rn^ 
is made equal to rriK- The required additive con- 
stants are thus determined. A completely reliable 
comparison between JkI f-w in full QCD and in 
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Figure 4. Monte Carlo data for the shift in beta 
of the deconfining transition from the valence 
approximation to full QCD, for Kogut-Susskind 
fermions, in comparison to the one loop weak- 
coupling prediction. 



the valence approximation depends on unknown 
higher terms in the chiral lagrangians in both 
cases. With plausible guesses for these terms, 
Sharpe [20] obtains the result that this quantity 
in the valence approximation will be smaller than 
in full QCD by about 0.12. The calculation of 
Ref. [24] finds Jk I f-w below its experimentally ob- 
served value by 0.06(9). 

Still another method of calculating differences 
between the valence approximation and full QCD 
has been considered by Hasenfratz and De- 
Grand [21]. These authors evaluate the effective 
shift in chromoelectric charge arising from vac- 
uum polarization due to heavy Kogut-Susskind 
quark-antiquark pairs in full QCD. For suffi- 
ciently heavy quarks, the determination of this 
change is reduced to evaluating the one loop 
quark-antiquark vacuum polarization term cou- 
pling two gluons in the weak-coupling expansion 
for lattice QCD. The resulting shift in (3 is linear 
in the number of quark flavors. The predicted A/3 
is compared with Monte Carlo data for the dif- 
ference between the critical (3c of the deconfining 



transition for QCD without quarks and (3c with a 
number of quarks ranging from 2 to 24. For lat- 
tices with time direction period ranging from 4 to 
8, Figure 4 shows the predicted shift in compari- 
son to Monte Carlo data. The overall agreement 
is quite good for TOjO > 0.05. 

2.4. Wave Functions, State Vectors 

I would like to mention briefly two recent pieces 
of work on hadron wave functions and state vec- 
tors. Kieu and Negele [22] have examined a 
Bethe-Salpeter wave function for the pion with 
quark and antiquark joined, in effect, by the 
ground state chromoelectric field configuration 
which would occur between a static color charge 
at the quark position and static anticolor charge 
at the antiquark position. They found that this 
wave function falls off with the quark-antiquark 
separation significantly more slowly than does 
the Coulomb gauge wave function. Thus for 
increasing quark-antiquark separation, a quark- 
antiquark pair joined by the ground state field 
configuration becomes progressively closer to the 
true pion state vector than does a Coulomb gauge 
quark-antiquark pair. 

Liu [23] has compared hadron masses calcu- 
lated in the usual valence approximation with 
masses obtained in the valence approximation us- 
ing a quark coupling matrix that allows quarks 
to propagate in only one direction of time. The 
state vector of a hadron in the usual valence ap- 
proximation, Liu shows, at any instant may in- 
clude quark-antiquark pairs in addition to the va- 
lence quarks. These pairs do not occur in hadrons 
composed of quarks which can propagate in only 
one time direction. As consequence of eliminat- 
ing these pairs, Liu finds, the rho and pion be- 
come nearly degenerate and the nucleon and delta 
baryon become nearly degenerate. 

3. Spectrum Calculations with Fixed Lat- 
tice Spacing 

3.1. Valence Approximation 

Three recent spectrum calculations using the 
valence approximation are summarized in Ta- 
ble 4. Kim and Sinclair used Kogut-Susskind 
fermions with wall sources and point sinks. The 




□ degenerate quarks 

X non— degenerate quarks 



0.25 

3.49 3.50 3.51 3.52 3.53 3.54 

l/(2«:.„) 



Figure 5. UKQCD data for vector mesons masses. 
For quark and antiquark hopping constants fei 
and ^2, Kh is (2fei)"^ + (2^2)"^ 



UKQCD collaboration used the clover action. 
The hadron propagators were found with point 
sources and point sinks. The clover action [29] 
consists of the Wilson fermion action with the 
simplest additional term added to cancel, to ze- 
roth order in the chromoelectric coupling con- 
stant, the 0(a) irrelevant coupling of the Wilson 
action. 

In addition to calculating pseudoscalar and 
vector masses with quark and antiquark masses 
equal, the UKQCD collaboration also found 
masses for mesons composed of a quark and an- 
tiquark with different masses. With this data a 
direct test can be made of Eq. (2) and the corre- 
sponding equation for the squares of pseudoscalar 
meson masses. Figure 5 shows a comparison of 
the masses of vector mesons composed of quark 
and antiquark with the same mass and the masses 
of vector mesons composed of quark and anti- 
quark with different masses. The data strongly 
supports Eq. (2). Corresponding data for pseu- 
doscalar masses supports the version of Eq. (2) 
for squared pseudoscalar masses. 





action 


lattice 


/3 


rUqa 


m^/mp 


configs. 


ref. 


Argonne 


valence, KS 


323 X 64 


6.0 


0.0100 


0.500(7) 


66 


[25] 






323 X 64 


6.0 


0.0050 


0.393(10) 


66 


[25 






323 X 64 


6.0 


0.0025 


0.311(9) 


66 


[25 






323 X 64 


6.5 


0.0100 


0.659(8) 


100 


25 






323 X 64 


6.5 


0.0050 


0.520(12) 


100 


[25 






323 X 64 


6.5 


0.0025 


0.401(15) 


100 


25 




action 


lattice 


/3 


k 


m^/rrip 


configs. 


ref. 


UKQCD 


valence, clover 


243 X 48 


6.2 


0.14144 


0.77(1) 


60 


[261 






243 X 48 


6.2 


0.14226 


0.62(2) 


60 


[261 






243 X 48 


6.2 


0.14262 


0.52(2) 


60 


[261 


GFll 


valence, Wilson 


163 X 32 


5.70 


0.1600 


0.856(3) 


219 


[61 










0.1650 


0.690(6) 


219 


[61 










0.16625 


0.612(6) 


219 


[61 










0.1675 


0.491(8) 


219 


[61 






243 X 32 


5.70 


0.1600 


0.854(2) 


92 


[61 










0.1650 


0.693(4) 


92 


[61 










0.1663 


0.610(5) 


58 


[61 










0.1675 


0.502(6) 


92 


[61 






243 X 36 


5.93 


0.1543 


0.830(5) 


210 


[61 










0.1560 


0.737(6) 


210 


[61 










0.1573 


0.603(6) 


210 


[61 










0.1581 


0.466(8) 


210 


[61 






322 X 30 X 40 


6.17 


0.1500 


0.867(2) 


219 


[61 










0.1519 


0.735(4) 


219 


[61 










0.1526 


0.633(6) 


219 


[61 










0.1532 


0.478(1) 


219 


[61 




action 


lattice 


/3 


k 


m^/rrip 


time 


ref. 


HEMCGC 


full QCD, Wilson 


163 X 32 


5.3 


0.1670 


0.722(4) 


2400 


[271 






163 X 32 


5.3 


0.1675 


0.599(8) 


1300 


[271 



Table 4 

Parameters of recent hadron spectrum calculations. 
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3.2. Full QCD 

Data for full QCD with two flavors of Wil- 
son fermions collected by the HEMCGC collab- 
oration [27] is also listed in Table 4. These 
calculations use the hybrid Monte Carlo algo- 
rithm. Hadron propagators were calculated for 
wall sources and either wall or point sinks, and in 
addition for Gaussian sources with either Gauss- 
ian or point sinks. The integrated autocorrelation 
time for the plaquette was found to be about 80 
time units for the data at the heavier pion mass 
and 120 time units for the data at the lighter pion 
mass. An examination of the dispersion of pion 
effective masses calculated from propagators av- 
eraged over bins gave roughly equivalent autocor- 
relation estimates. Thus the data collected at the 
larger pion mass consisted of perhaps 30 indepen- 
dent configurations and the data with the smaller 
pion mass perhaps 10. Figures 6 and 7 show full 
QCD results obtained with various combinations 
of sea quark and valence quark hopping constants 
in comparison to valence approximation masses 
for Wilson fermions on a lattice 16^ x 32 at /3 of 
5.85 and 5.95 [28]. For the range of parameters 
considered, the differences between the valence 
approximation and full QCD consistent with the 
statistical uncertainties. 



2.0 r^ 
1.8 \- 

g°- 1.6 \- 

B 1.4 - 
1.2 - 



1.0 



0.0 



O.D 
m„/m„ 



Figure 6. HEMCGC mass ratios with various 
combinations of valence quark and sea quark hop- 
ping constants. The squares have sea quark k of 
0.1675, the crosses have sea quark k of 0.1670. 
The diamonds are valence approximation results. 
The circle and question mark show expected val- 
ues for infinite quark mass and for the real world 
value of the light quark masses. 



4. Valence Approximation Zero Lattice 
Spacing Limits 

4.1. Masses 

The first systematic evaluation of the infinite 
volume continuum limit of hadron masses ex- 
trapolated to physical quark masses was reported 
this year by the GFll collaboration [6]. Hadron 
masses were calculated with Wilson quarks on a 
set of four lattices. The parameters entering these 
calculations are listed in Table 4. 

Each configuration was fixed into Coulomb 
gauge. In all cases, hadron propagators were then 
calculated using a Gaussian quark source with 
mean radius squared < r^ > of 6 in lattice units, 
and using point sinks and Gaussian sinks with 
< r^ > of 1.5, 6, 13.5, and 24. Hadron masses 
were fitted to each channel by examining effective 
mass plots and selecting the largest interval which 
might appear to be a plateau. An automatic fit- 




0.1 0.2 0.3 

(m -m„)/(3m + m„) 



Figure 7. HEMCGC hyperfine splittings. Sym- 
bols have same meaning as in preceding figure. 
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ting program then did a fit to each subwindow of 
four or more time slices within this region. The 
subwindow giving the smallest value of y} per de- 
gree of freedom was selected as the final fit. Fits 
to each window were done by minimizing the fit's 
correlated y} among all data points fitted. Fits 
were also done in this way simultaneously to the 
propagators for point sinks and Gaussian sinks 
with < r^ > of 1.5 and 6. Error bars on all of 
these results were found by the bootstrap method. 
A examination of effective mass plots showed, 
in all cases, that the propagators for hadron sinks 
with larger < r^ > reached plateaus at smaller 
time separations between source and sink and 
therefore coupled more weakly to excited states. 
Thus a comparison of masses fitted to channels 
with different radii provided a test for the pres- 
ence of high mass contamination. Typical ef- 
fective mass plots and fits are shown in Fig- 
ures 8, 9 and 10. These figures show the rho 
on 30 X 32^ x 40 with k of 0.1532, for sinks with 

< r^ > of 0, 6 and 24. Table 5 gives the masses 
obtained from these fits, along with masses from 
fits to sinks with < r^ > of 1.5 and 13.5, and 
the mass from a simultaneous fit to data with 

< r^ > of (point sinks), 1.5 and 6. The con- 
sistency of these results is evidence that the low- 
est mass has been picked out accurately in each 
channel. Similar results are obtained from prop- 
agators found on the other GFll lattices listed in 
Table 4. Thus the masses determined from these 
fits appear to be free of the high mass contami- 
nation which Ukawa [5] suggested might occur. 

For each GFll data set of Table 4, the critical 
hopping constant fee was determined, as usual, by 
fitting m\ in lattice units to C[(2fe)~-^ - (2fec)~'^]- 
With quark mass TOj then defined as {2\i)~^ — 
{2k ^~^, hadron masses for a continuous range 
of quark masses were found by fitting the Monte 
Carlo data to linear functions of quark mass. 
Figure 11 shows such fits to hadron masses at 
the three lightest values of TOj for the lattice 
30 X 32^ X 40 at /3 of 6.17. The masses were de- 
termined by simultaneous fits to propagators with 
sink < r^ > of 0, 1.5 and 6. The up and down 
quark masses were then taken to be equal and 




Figure 8. Effective masses and fitted mass for 
the rho propagator with point sink on the lattice 
30x32^x40 with k of 0.1532. Dotted vertical lines 
bound the window selected by the automated fit- 
ting procedure. 
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given by the "normal" quark mass. 



which 



Figure 9. Effective masses and fitted mass for 
the rho propagator with Gaussian sink < r^ > of 
6. on the lattice 30 x 32^ x 40 with k of 0.1532. 
Dotted vertical lines bound the window selected 
by the automated fitting procedure. 



was determined by fitting the ratio of extrapo- 
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Figure 10. Effective masses and fitted mass for 
the rho propagator with Gaussian sink < r^ > of 
24. on the lattice 30 x 32^ x 40 with k of 0.1532. 
Dotted vertical lines bound the window selected 
by the automated fitting procedure. 



hadron 


<r^ > 


mass 


range 


x' 


pion 





0.1445(20) 


12-15 


0.2 




1.5 


0.1453(20) 


15-18 


0.4 




6 


0.1451(19) 


13-17 


0.1 




13.5 


0.1453(20) 


13-16 


0.0 




24 


0.1455(20) 


13-16 


0.0 


rho 


0,1.5,6 


0.3024(48) 


13-16 


0.4 







0.2991(59) 


13-16 


0.2 




1.5 


0.2987(59) 


12-15 


0.1 




6 


0.2969(66) 


13-16 


0.1 




13.5 


0.2965(76) 


13-16 


0.0 




24 


0.2965(88) 


13-16 


0.0 


nucleon 


0,1.5,6 


0.4234(65) 


13-16 


1.1 







0.4187(59) 


12-15 


0.5 




1.5 


0.4188(87) 


12-17 


0.7 




6 


0.4199(88) 


13-16 


0.5 




13.5 


0.4215(80) 


13-16 


0.2 




24 


0.4097(78) 


10-13 


0.9 


delta 


0,1.5,6 


0.5267(81) 


13-17 


0.9 







0.5233(88) 


14-17 


0.1 




1.5 


0.5188(111) 


14-17 


0.0 




6 


0.5133(104) 


14-17 


0.5 




13.5 


0.5197(101) 


11-16 


0.9 




24 


0.5120(101) 


11-15 


0.7 


Table 5 










Fitted masses from the lattice 32 
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with k of 0.1532 
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Figure 11. For a 30 x 32^ x 40 lattice at (3 of 6.17, 
TO^, rrip, rriN and toa, in units of the physical 
rho mass TOp(to„), as functions of the quark mass 
rUq , in units of the strange quark mass rUs . The 
symbol at each point is larger than the error bars. 



lated values, rriT^/mp to experiment. The strange 
quark mass was found using the pseudoscalar ver- 
sion of Eq. (2) to determine rriK, then fitting 
the ratio of extrapolated values mK/mp to ex- 
periment. Masses for other strange hadrons were 
found using Eq. (2) and its analogues for baryons. 

For Wilson fermions, hadron mass ratios are 
expected to approach their continuum limits as 
linear functions of lattice spacing. Linear extrap- 
olations of hadron mass ratios to zero lattice spac- 
ing are shown in Figures 12 - 14. The masses in 
this extrapolations were found from simultaneous 
fits to sinks with < r^ > of 0, 1.5 and 6 for the 
runs on 16^ x 32 with j3 of 5.7, 24^ x 36 with j3 of 
5.93 and 30 x 32^ x 40 with beta of 6.17. The val- 
ues of /3 in each case are chosen so that the lattice 
period L is the same in physical units in all cases 
and within statistical errors equal to 9to~-'-. 

The continuum ratios found in finite volume 
were then corrected to infinite volume. After 
extrapolation to physical quark mass, hadron 
masses measured in units of rUp on the lattice 
16^ X 32 were compared with ratios on the lattice 
24^ X 32. The changes were all found to be less 
than 4.5%. Using either Luscher's formula [11] 
for the approach of masses to their infinite vol- 
ume limit or the model of Ref. [10], it can then 
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1.40 




Figure 12. Mass ratios and linear fits as functions 
of the lattice spacing a in units of l/rUp. The 
error bars at zero lattice spacing are uncertainties 
in the extrapolated ratios, and the points at zero 
lattice spacing are observed values. 
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Figure 13. Mass ratios and linear fits as functions 
of the lattice spacing a in units of l/rUp. The 
error bars at zero lattice spacing are uncertainties 
in the extrapolated ratios, and the points at zero 
lattice spacing are observed values. 




Figure 14. Mass ratios and linear fits as functions 
of the lattice spacing a in units of l/rUp. The 
error bars at zero lattice spacing are uncertainties 
in the extrapolated ratios, and the points at zero 
lattice spacing are observed values. 

be argued that masses on 24^ are within 1% of 
their infinite volume limit. 

For any hadron mass m measured in units of 
rUp define the finite volume correction term for 
lattice period L 



A{(3,L) = —{(3,<^) 
m„ 



<I3,L). 



(8) 



The quantity needed is A(oo, Qto"-*-). Based on 
the size of the change in m/mp from /3 of 5.7 to 
infinity and on the size of the change in this ratio 
from 24^ x 32 to infinite volume, it is then argued 



A(oo,9m;i)«— (5.7,9m;i)- 
mp 

— (5.7,13.5m7i), 



(9) 



to within an error of less than 2% of m/mp. In- 
finite volume continuum mass ratios were then 
found using Eq. (9). 

Table 6 lists continuum hadron mass ratios 
found in finite volume and after applying the 
infinite volume correction. The infinite volume 
numbers taken together are statistically consis- 
tent with experiment. The differences between 
these numbers and experiment are less than 6%, 



14 



ratio 


fin. vol. 


inf. vol. 


obs. 


rriK'/mp 


1.149(9) 


1.166(16) 


1.164 


m^/mp 


1.297(19) 


1.333(32) 


1.327 


rriM/mp 


1.278(68) 


1.216(104) 


1.222 


Am/mp 


1.865(45) 


1.927(74) 


2.047 


rriA/mp 


1.609(78) 


1.565(122) 


1.604 


rriY;' Imp 


1.805(58) 


1.806(80) 


1.803 


m^. /nip 


2.009(49) 


2.055(65) 


1.996 


mn/mp 


2.205(56) 


2.296(89) 


2.177 


AS^/™. 


0.305(8) 


0.319(12) 





decay fin. vol. inf. vol. 



obs. 



Table 6 

Calculated hadron mass ratios in comparison to 

observed. The mass difference Arn is toh + tos — 

mjv. 

equivalent to less than 1.6 times the statistical 
error. 

The value of A^— — /toq shown is obtained 

MS' ^ 

by extrapolation to zero lattice spacing of 
(^hX-L-a)/{mpa) as a linear function of rUpa. Here 

h!^-L-a is determined from a-rr^ by the two-loop 

MS . . MS J ^ ^ 

Callan-Symanzik equation, and a-Ms ^® found 
from /3 following Ref. [30]. The slope of the fit 
of (A^— ^a)/(TOpa) to a linear function of rUpa is 
statistically consistent with zero. Thus iripa fol- 
lows asymptotic scaling within statistical errors. 

4.2. Decay Constants 

The first evaluation of the infinite volume con- 
tinuum limits of meson decay constants was re- 
ported by the GFll group in Ref. [24]. The cal- 
culation used the same the data sets listed in Ta- 
ble 4. The fits to propagators and extrapolations 
to physical quark mass, zero lattice spacing and 
infinite volume were all done by versions of the 
methods used in the mass evaluation of Ref. [6]. 
Table 7 shows predicted results in comparison to 
experiment. The finite renormalizations factors 
in these numbers are found to one-loop order fol- 
lowing the mean-field improved perturbation ex- 
pansion of Lepage and Mackenzie [31]. Statis- 
tically equivalent results but with slightly larger 
error bars are given by zero order mean-field im- 
proved finite renormalizations. Equivalent results 
with still larger error bars are given by naive finite 
renormalizations. 

Although overall the predicted values are in 
fair agreement with experiment, the circumstance 



U/rUp 0.106(9) 

fK/rUp 0.121(6) 

Fp/mp 0.177(21) 

F^/mp 0.217(19) 



0.106(14) 0.121 

0.123(9) 0.148 

0.173(29) 0.199 

0.253(35) 0.219 



Table 7 

Calculated values of meson decay constants ex- 
trapolated to zero lattice spacing in finite volume, 
then corrected to infinite volume, compared with 
observed values. 

that three of four decay constants range from 
0.9 to 2.8 standard deviations below experiment, 
suggests that valence approximation decay con- 
stants may lie systematically somewhat below ex- 
periment. A simple physical argument supports 
this possibility [30]. The valence approximation 
can be thought of as the replacement of the mo- 
mentum and frequency dependent color dielectric 
constant caused by quark-antiquark vacuum po- 
larization with its low momentum limit [13]. At 
low momentum the effective quark charge in the 
valence approximation agrees with the low mo- 
mentum effective charge of the full theory. The 
valence approximation could therefore be fairly 
reliable for low-lying baryon and meson masses 
if, as seems likely, these masses are determined 
largely by the long distance behavior of the chro- 
moelectric field. Due to the absence of dynamical 
quark-antiquark vacuum polarization, however, 
the quark charge in the valence approximation 
falls faster with momentum than it does in the 
full theory. As a result at short distance the at- 
tractive quark-antiquark potential in the valence 
approximation is weaker than in the full theory, 
and meson wave functions in the valence approx- 
imation are smaller at the origin less than in the 
full theory. Meson decay constants are propor- 
tional to the square of wave functions at the origin 
and thus should be to be smaller in the valence 
approximation than in the full theory. 

4.3. Glueball Masses 

An evaluation of the infinite volume, contin- 
uum limit of glueball masses was reported by a 
Liverpool- Wuppertal collaboration [32]. The lat- 
tice sizes, ensemble sizes and (3 of the new calcu- 
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lattice 


/3 


configs. 


TOo++a 


Liverpool- Wuppertal 


16^ X 32 
32^ 
32^ 


6.0 
6.2 
6.4 


3040 
1001 
3200 


0.701(27) 
0.522(38) 
0.415(14) 


GFll 



16^ X 24 


5.70 


8094 


0.983(40) 


20^ X 30 


5.83 


4002 


0.858(43) 


16^ X 24 


5.93 


30640 


0.786(12) 


123 X 24 


5.93 


48278 


0.771(10) 


24^ X 36 


6.17 


31150 


0.582(10) 


30 X 322 X 40 


6.40 


25440 


0.433(11) 



Table 8 

Parameters of recent glueball spectrum calcula- 
tions. 

lation are shown in Table 8. This calculation used 
the same methods as earlier work by Michael and 
Teper at smaller values of /3 [33]. Glueball oper- 
ators were constructed by an iterative gauge in- 
variant smearing process which, at each step, re- 
turns smeared links which have twice the length 
of the links taken as input [34]. For each time 
separation between source and sink operators, a 
matrix of correlation functions was found among 
the various operators carrying the same quan- 
tum numbers. The possible operators for each 
set of quantum numbers were made by summing 
together rotations of loops with the same shape. 
The maximal eigenvalue of this matrix was then 
extracted. The logarithm of the ratio of maximal 
eigenvalue at separation t to the maximal eigen- 
value at separation t + 1 defines the eff'ective mass 
at separation t. By including a range of different 
shapes of loops in the set of operators, then max- 
imizing the eigenvalue of the correlation matrix, 
combined operators were produced which appear 
to couple very efficiently to the lightest state in 
each channel. An comparison of effective masses 
produced this way and with simpler choices of 
glueball operator suggests that simpler choices do 
not select out the lightest states as efficiently. 

The final mass value in each channel was found 
generally by looking for the smallest time sepa- 
ration giving an effective mass statistically con- 
sistent with the effective masses at larger time 
separation. Figure 15 shows effective masses and 




Figure 15. Effective masses and final mass value 
found by the Liverpool- Wuppertal group for the 
0++ glueball on a lattice 32^ at j3 of 6.4. 

the mass taken as the final value for the 0"*""*" glue- 
ball at /3 of 6.4. In this case, the effective mass 
between t of 2 and 3 was chosen as the final mass 
value. Although this choice appears reasonable, it 
is not clear whether the data has actually reached 
a plateau at this point. It appears possible that 
the rate of growth in error bars may be obscuring 
a continued decline in the effective mass. A large 
collection of masses were also calculated for other 
glueball states. The effective mass plots again did 
not show clear plateaus. The final masses in most 
other channels were therefore taken only as upper 
bounds. 

To evaluate the continuum limit, mQ++a in lat- 
tice units was divided by the square root of the 
string tension in lattice units \fKa. This ratio 
is expected to approach its continuum limit lin- 
early in a^. A linear fit was made to the ratios 
(TOo++a)/(vira) as a function of Kc^ . The con- 
version of mass ratios to MeV was done using 
an estimate of K in physical units of 440 MeV, 
based on a variety of indirect arguments. The 
value of TOo++ at /3 of 6.4 was found to be about 
50 MeV below the continuum limit of 1600 MeV. 
The mass at 6.4, 1550 ±50 MeV, was then quoted 
as the continuum limit with a additional system- 
atic extrapolation error of 50 MeV. No estimate 
was given for the error associated with the choice 
of value of K. 

An alternate way to obtain a continuum 
limit prediction for TOo++ is to fit the ratio 
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Figure 16. Liverpool- Wuppertal data for 

(TOo++a)/(A^— ^a) and a linear extrapolation to 
zero lattice spacing. 

(TOo++a)/(A^-^a) to a linear function of (A^-^a)^. 

A value of A^-^ in physical units is given by the 

continuum limit of hX-l—/mp reported in Ref. [6]. 
This value is free of the ambiguities in the estima- 
tion of the physical value of K. Figure 16 shows 
the continuum limit of TOo++ obtained in this way. 
Data at lower (3 is also included from Refs. [33] 
and [35]. The predicted continuum too++/A^-^ is 

6.67±0.33. The value of A^-^ from the continuum 

MS 

limit of {A^°La)l{m^a) is 243.7±6.8 MeV. This is 
statistically consistent with the value found using 
TOp but has a somewhat smaller error bar. Includ- 
ing the uncertainty in A^— ^, the final prediction 
obtained from the Liverpool- Wuppertal data for 
mo++ is then 1625 ± 92 MeV [36]. 

A calculation of the infinite volume continuum 
limit of glueball masses was also reported by the 
GFll collaboration [37]. The lattice sizes and (3 
used are shown in Table 8. For the lattice 20^ x 30 
results are shown for smeared glueball operators 
defined in Coulomb gauge. For the four other 
lattices gauge invariant smeared glueball opera- 
tors are constructed following Ref. [38]. Smeared 
links are defined by adding to each link the sum of 
space direction staples with a coefficient e which 
was taken to be 0.25 at (3 of 5.7 and 1.0 other- 
wise. This step is then iterated some number of 
times Ns- From the smeared links, square loops 
Nl X Nl were constructed. Sums over rotations of 
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Figure 17. Effective masses and final mass value 
found by the GFll group for the 0"*""*" glueball on 
a lattice 30 x 32^ x 40 at /3 of 6.4. 



such loops then give operators for the 0"*""*", 1"' 
and 2+"'" glueball. A correlation matrix among 
the operators with each set of quantum num- 
bers for a range of Ns and Ni, was evaluated at 
each time separation. Effective masses were de- 
termined from the diagonal terms in these matri- 
ces and intervals of smearing sizes and loop sizes 
found for which effective mass plots showed clear 
plateaus. Figure 17 shows the effective mass plot 
for the 0++ glueball at /3 of 6.4 found for Ns of 8 
and Nl of 9. A clearer plateau occurs here than in 
Figure 15 since the ensemble used for Figure 17 is 
about eight times larger. Similar plateaus in 0"*""*" 
effective masses were found for several other com- 
binations of Ns and Ni, for this value of (3 and 
for the other values listed in Table 8. A mass was 
extracted from the propagator of Figure 17 fitting 
the time interval indicated. Masses were found in 
other channels similarly. An optimal mass was 
then found among the set of channels showing 
plateaus using a correlation matrix among these 
masses generated by the bootstrap method. Al- 
ternative fits using the full underlying propagator 
matrix gave statistically equivalent results. 

For the 2+"'" the effective mass plateaus were 
weaker than those for the 0"*""*" and below (3 of 
5.93 a clear signal could not be found. For the 
1"' glueball an acceptable signal could not be 
found in any of the data sets of Table 8. 

Figure 18 shows the masses found for 
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Figure 18. GFll date for (mo++a)/(A^a) and 
a linear extrapolation to zero lattice spacing. 
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7.14 ± 



(TOo++a)/(A^— — a) to zero lattice spacing. 

predicted continuum too++/A^— — is 

0.21. Using A^— — from the continuum limit of 

^ MS 

(A^— — a)/(TO^a) and including the uncertainty in 

A^, the final predicted mo++ is 1740 ±71 MeV. 

A comparison of the GFll data for the lat- 
tice 16^ X 24 and 12^ x 24 at /3 of 5.93 com- 
bined with Luscher's formula for volume depen- 
dence [11] gives a bound of 0.5% on the difference 
between this prediction and its infinite volume 
limit. A similar estimate applies to the Liverpool- 
Wuppertal value of TOo++ . 

The Liverpool- Wuppertal values for TOo++ at 
fixed lattice spacing, and for the continuum limit 
both in units of A^— — and in MeV are statisti- 

MS 

cally consistent with the GFll values. The GFll 
values, however, appear to be more reliable be- 
cause of the larger ensembles used and the result- 
ing more clearly defined effective mass plateaus. 
The two well established 0"*""*" states within any 
reasonable distance of 1740±71 MeV are /o(1590) 
and /o(1710). The first of these is 2 standard de- 
viations low thus /o(1710) appears to be favored. 
A question which remains is why a valence ap- 
proximation glueball calculation should be be- 
lieved at all. A commonly held opinion is that 
the 0"*""*" may mix strongly with quark-antiquark 
states and might be shifted by as much as 20% in 



going to full QCD. My guess, on the other hand, 
is that the mass shift due to mixing is likely to 
be less than 4%. A 0"*""*" quark-antiquark state 
must have odd orbital angular momentum. The 
resulting angular momentum barier should tend 
to suppress the quark-antiquark annihilation re- 
quired for mixing with chromoelectric field. An 
estimate of the effect on TOo++ of mixing with 
quark-antiquark states is given by the mass shift 
in quark-antiquark states, with nonzero orbital 
angular momentum, due to mixing with chromo- 
electric field. In the valence approximation the 
isosinglet 2+"'" state /2(1270) and isovector 2+"'" 
state 02(1320) are exactly degenerate. In the real 
world, the /2 mixes with glue and the 02 does 
not. Thus the splitting between the /2 and 02 is 
the shift in the /2 due to mixing. This splitting, 
and therefore the shift, is 43 ± 5 MeV. Similar 
estimates can be made comparing /ii(1170) with 
61(1235), /i(1285) with ai(1260), W3(1670) with 
P3(1690), and perhaps also comparing /o(975) 
with ao(980) if these are quark-antiquark states 
rather than kaon molecules. The splittings be- 
tween these pairs range from 9±3 MeV to 60 ±20 
MeV. Thus it seems to me reasonable to guess 
that the shift in the mass of the 0"*""*" glueball due 
to mixing with quark-antiquark states will be less 
than about 60 MeV, which is 3.4%. 

The other possible source of error in the va- 
lence approximation to the 0"*""*" mass arises from 
virtual quark-antiquark loops contributing more 
than just a shift in the effective chromoelectric 
charge. This same error, however, is present 
in the valence approximation to the masses of 
hadrons composed of quarks and antiquarks. The 
effect of this error on low-lying hadron masses, 
discussed in Sect. 4.1, is less than 6%. This error 
combined with the effect of mixing could decrease 
1740 ±71 MeV enough to make it consistent with 
/o(1590). The /o(1710) would still appear to be 
favored. 

I am grateful to a large number of people for 
answering my email, explaining their work, and 
sending data and pictures. I would particularly 
like to thank G. Bali, T. DeGrand, S. Gottlieb, 
Y. Iwasaki, T. Kieu, Y. Kim, J. Labrenz, D. 
Leinweber, K.F. Liu, J. Negele, F. Rapuano, D. 
Richards, K. Schilling, S. Sharpe, A. Ukawa, A. 
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